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ABSTRACT 


The purpose of this project is to localize buried 
objects in water-saturated sand by variable incidence 
acoustic pulse reflections. In particular, thin aluminum 
plates are used to model the buried objects. The sound 
speed and the depth of the overlying sand layer are predicted 
by using the Wee ee and the ray parameter methods. The 
experimentally computed results have an average 3.5% depth 
error. The calculated speed of sound in the water-saturated 
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sand layer differed by 8.1% between the T?-x and the ray 


parameter methods. 
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I. INTRODUCTION 


A. OBJECTIVES 

Identification of buried objects in the ocean bottom is 
an important problem in many fields such as oil exploration, 
geophysical research, and ocean engineering. A simpler 
problem is the identification of layers of infinite extent 
and unknown properties in the ocean sediments. The objec- 
tives of this thesis are: 


1. to determine experimentally the depth of burial 
of a thin plate in a simulated sedimentary bottom; 


2. to determine the sound speed and the thickness 
of the overlying sand layer; and, 


3. to determine the sound speed and thickness of 
the thin plate. 


It was the intent of this thesis to determine all of these 
properties by remote sensing with acoustic pulses in order 
to simulate an ocean measurement with a ship and near- 
surface sensors. 

The first objective could have been accomplished by 
normal-incidence time-of-flight measurements; however, the 
second objective requires more information than that so a 
variable angle of incidence technique was used. This second 
objective was actually the principal goal since the third 
objective could not be achieved by using the technique 
described herein. The technique failed because the pulse 


could not be made short enough to distinguish the reflection 


from the plate’s lower surface. This is primarily the 

result of using very thin plates with respect to the acoustic 
wavelength. Perhaps the amplitude of the reflection from 

the top of the plate could be used to assist in idenerives 
the plate material but the overlying sand layer makes this 


VEry acter iculc. 


B. APPROACH 


The approach used is as follows: 


1. Short pulses consisting of several cycles of a 
150kHz sinusoidal wave were used as the source 
pulse. 

2. A directional projector was used to minimize inter- 


ference from surface and side-wall reflections. 


3. The receiving hydrophone was mounted on a carriage 
so that the horizontal distance from source £6 
receiver could be easily varied. 


4. Each of the thin plates were buried in the sand 
bottom of a small water tank. 


5. The plates were aluminum ranging in thickness from 
0.1 to 1.8 cm and were buried at depths from 1.8 
to Sat@em. 

6. Two analytical methods were used for reducing the 


data. One method was based on the relationship 
between the square of the travel time (T) of the 
reflections and the square of the source-to-receiver 
horizontal range (X). The other method was based 
on maintaining a constant ray parameter (P) defined 
as the sine of the incident angle divided by the 
local sound speed. 
C. BACKGROUND 
Reference 1 gives an excellent historical development 
of various seismological methods for determining properties 


of sediment layers. The following is a brief summary of 


that eccoune. 
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Karcher and HaSeman made the first successful seismic 
reflection measurements in 1921. Techniques for acoustic 
determination of sediment properties have advanced greatly 
Since then and generally fall into one of three categories: 
normal incidence, wide-angle reflection, and refraction 
methods. Most of the development has been done by the 
petroleum exploration industry; however, the geophysical 
community has used these methods to study the earth's crust 
as early as the Tate 1930's. Dix, Krey and Durbaum 
refined the sound speed measurements for sediment layers 
below the first layer by properly accounting for the shape 
of the arrival time as a function of range curves. By 
1955, Dix had demonstrated the equivalence of T-X methc 


Mei <= <* meted and, in’1968, Le Pichon refined the 


~?-x? method until it became a standard analysis tool. 
A major enhancement was introduced by Bryan (see Ref. 2) in 
1974 with the ray parameter method. Both the ee method 


and the ray parameter method were used in the data analysis 


for this thesis. 


ie 


ti. THEGRS 


A. T=-X AND T?-x? METHODS 


A simple model of a layered sediment is shown in Fig. 
2.1. The layers have thicknesses Has H. and H3 and sound 
1’ V5 and V3- The T?-x2 method [Ref. 1] can be 
applied to this situation as described in this section. 


speeds V 





WATER 
/ 
\ 
SAND Hy. 
f 
0,3 PLATE : 
Gen 
Figure 2.1. Geometry of Layers 


For convenience, let the sound speed in the uppermost 


layer (the water layer) be less than that in the first 


le. 


sediment layer. The length of the path that reflects from 


the top of the sediment layer is given by, 


2 a Z ee 
Rui —) >.0m 5 4(H, Ss) (1) 


To convert this to travel time squared, divide through 


by Vi: 


T*(X) = T°(0) + (—5)x* (2) 
a 

where T(X) is the reflection travel time at horizontal range 
X and T(0) is the intercept of the travel time curve at 
X = 0. Since T* (0) is a constant, the reflection curve is 
a hyperbola. The corresponding curve for the direct path 
from the source to the receiver is a straight line and this 
hyperbola is asymptotic to that direct-path line. 

When plotted in the Tie plane, however, these lines 
are both straight. In this way, the slope of the reflection 
curves can be easily determined. As can be seen from the 

15 eye 


equation for the reflection path, the slope in the T’-xX 


plane is + which, in this thesis work, is calculated by 


5 
pal 2.2 


linear regression of the T’-X” data. This rather simple 
but accurate method is due geneen [Ref. 1] and is asso- 
Ciated with the wide-angle seismic reflection technique. 

The paths that pass through one or more sediment layers 
before reflection (see, for example, Ri>5 Mee 2 ea Not 


lead exactly to hyperbolas but the departure is not great 
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Figure 2.2. T=<sand [exe Curves 
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and only a small amount of error is introduced by fitting 


2_x? plane. 


this data with straight lines in the T 

If the layers were extremely thick or had great differ- 
ences in sound speed, then the error might become apprecia- 
ble. The reciprocal of the slope of the R) 5 curve is equal 


to the average speed (Squared) through both the water and 


the first sediment layer: 


< 
| 


= (Vj T)+V5T,)/(T,+T,) (3) 


Representative T-X curves for this case are shown in 
Fig. 2.2.a. Curve D represents the direct path and Ri4 
represents the reflection curve for reflection from the top 


of the sediment. For R the exact curve (solid line) is 


daa 
given as well as the hyperbola approximation (dotted curve) 
to that curve. The equivalent curves in ene Space are 
Shown in Fig. 2.2.b. Here the inverse slopes give the sound 
speed squared. The intercept T(0) gives the two-way 
vertical travel time through the layers. 

The principal limitation of the Oa method is shown 
in these figures. For any reflection other than the reflec- 
tion from the top of the sediment, the curve is not 
eeeicily a straight line in noe coordinates but it falls 


Slightly below the equivalent straight line as range 


increases. 


arS 


B. THE RAY PARAMETER METHOD 

In order to measure the properties of a layer without 
knowing the structure of the overlying layers, the ray 
parameter method can be used [Ref. 2]. Any ray path is 
specified by a fixed value of the ray parameter P defined 


by, 
P = Sin €(2)7el2} (4) 


where z is depth, c(z) is the sound speed profile and @(z) 
is the angle that the ray makes with the vertical at depth 
Zz. Figure 2.3 shows the relationship between the incident 
angle and the ray parameter. 

According to the ray parameter method, the increment of 
travel time along the ray path gives the corresponding 


increments in depth and range: 
dz = cdtTcos 6 = (1l1-Pc) ord (5) 
qdX = cdtTcos @ = Pec'dt = Pc(l1l-P 


Pd (6) 


After integrating over depth from the surface to depth 


Zar these expressions can be rewritten as follows, 


az (7) 
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Figure 2.3. Geometry of Ray Parameter 


iese Cquations give the travel time T and the hori- 
zontal distance X along the ray path given by P in terms of 
Mmeeatals over depth. The equations must be solved for the 
sound speed variation c(z). This can be done for some 
Special cases of layer structure. 

Consider a homogeneous layer (sound speed cy) extending 
eEOm 2 


1 Boma wet Zeid. 2.5). hOr Ene ray Of a particular 


1 


ray parameter value, there are two (in this case) rays that 


ty, 


return from the bottom with that same parameter value 
although they arrive at the receiver at different ranges. 
Note that these rays are parallel through the water column. 
The travel time and horizontal distance for these two rays 


are, according €© Eqs= (7) andi 


1 7 —l 72 
AT(P) = 2 f en (1 =P7c7) eee (9) 
aa 
Z,+Az ne yD 
AXE ee Pc (1 -P“c*) dz (10) 
Z 
1 
where the factor 2 accounts for the round trip Of =he Zaye 
If the sound speed is constant at Cy in the layer, these 
equations can be integrated to give, 
7 =172 
AD = 2eea(1 gc) | Qe Gini 
5 2 m7 2 

AX = 2Pc, (1 -P C}) Az (125) 
The expressions for sound speed and layer thickness are 
then, 

u lee avi 
L72 
E> ft! Zee 
AZ = 7 C, AT(L -P Cy) (14) 
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Figure 2.4. T-X Curve and the Ray Parameter 


The important feature of the ray parameter method is that 
the structure above the layer has dropped out because of 
the choice of ray path. 

imorder t© apply this technique in practice, an initial 
value of ray parameter was chosen based on a geometry that 
gives good clean reflection signals. After the T-X curves 
are drawn, this value of ray parameter, which gives the 
Slope of the T-X curves at the points of interest, was used 
to identify the two parallel return rays. The principal 
difficulty of this method is that the T-X curves are not 


Straight lines and, therefore, curves must be fit to the 


ib) 


experimental data. This introduces some uncertainty into 
the process but, fortunately, the values for layer thickness 


is not very sensitive to the type of fit that is done. 
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III. EXPERIMENTAL PROCEDURES 


A. ARRANGEMENT OF THE EXPERIMENT 

The measurements were made in a tank which was 300 cm 
long, 115 cm wide and 90 cm deep. The tank was filled with 
carefully de-aerated (fresh) water and had a 55 cm deep 
sand bottom (also de-aerated). 

At first, an ITC 501 projector was used but the beam 
pattern was fairly broad leading to excessive reflections 
from the surface and the side-walls and, more importantly, 
this projector was large enough so that it would ring fora 
Significant length of time compared to the desired pulse 
length. As a result, the output pulse was too long regard- 
less of the length of the pulse used as input. 

These problems led to the selection of a USRF F-4l 
projector, which is a directional laboratory standard trans- 
ducer used in acoustic reciprocity calibrations over the 
frequency range from 15 to 150 kHz. This projector consists 
of lead zirconate-titanate elements cemented to a high- 
density kennametal disks. The resulting array of elements 
is about 3.8 cm wide and 5 cm high. 

The receiving hydrophone used was an omnidirectional 
LC-10 hydrophone mounted on a carriage that was free to roll 
along tracks installed on the long sides of the tank. In 


this way the horizontal distance (X) could be easily varied. 


Zedk 


A block diagram of the electronics used is shown in 
Fig. 3.1. The digital oscilloscope was used to measure the 
time differences between various pulses. In order to reject 
unwanted noise, a filter was used to pass only 100 kHz to 
200 kHz. 

The choice of pulse length was governed by the criteria 
that it be (1) long enough to reach steady-state conditions 
in the projector and (2) short enough to resolve the indi- 
vidual reflections from the layers. The best compromise 
turned out to be four cycles of a 150 kHz sine wave, but, 
unfortunately, this was not short enough to completely 
define the thin plate. At 150 kHz, one cycle is 6.7 uS 
so the four-cycle pulse was about 2/7 yS long. Amey eel 
received signal using the four-cycle pulse is shown in 
Fig. 3.2 while Fig. 3.3 shows the same situation with an 
eight-cycle pulse. The four-cycle pulse gives much better 
definition of individual arrivals. 

Before the actual measurements were made, several tests 
were performed in order to insure that side-wall and surface 
reflections would not interfere with the measurements. The 
depths of the projector and hydrophone were determined by 
trial and error not only to minimize interfering reflections 
but also to optimize the separation between the desired 
arrivals. In addition, this set a practical 1imicGyemeen 
horizontal range. 

An attempt was made to move the axis of the source-to- 


receiver line toward the tank wall so that two parallel 
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Figure 3.2. 4 Cycle Pulse Signal 





Figure 3.3. 8 Cyele Pulse "Siena 


measurements could be made without disturbing the sand to 


bury each plate. This introduced unacceptable side-wall 
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Figure 3.4. Location of Projector, Hydrophone 
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Figure 3.5. Projector Is Not on the Center Line 
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reflections so, instead, three plates were buried in line 
along the long axis of the tank with the projector and 
receiver positioned as shown in Fig. 3.4. 

Here, the plates were buried at different depths and the 
projector and hydrophone were positioned as shown. (Notice 
that in the third set-up on the right the source and receiver 
are reversed to avoid reflections from the end-wall.) 

Figure 3.5 shows the side-wall reflections that are 
introduced when the source-to-receiver axis is displaced 
only 15 cm toward one side of the tank. This result should 
be comparable to Fig. 3.2 but only a position along the 


centerline of the tank turned out to be acceptable. 


B. EXPERIMENTAL PROCEDURE 

Once the experiment was set up, the first step was to 
identify the individual arrivals from a particular buried 
plate. By considering the projector depth S, the hydrophone 
depth h, the water depth H, the major incident angle, and 
the slant distance R, the time of arrival of the various 
reflections was determined by simple image theory [Ref. 3, 
pp. 427-429]. A representative diagram of the various 
arrivals is shown in Fig. 3.6 for constant sound speed in 
the water layer. 

Since a directional projector was used and was aimed 
downward toward the sediment, the surface-reflected path 
7 could be neglected. For most of the measurements, paths l, 


3 and 4 were used. 
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Figure 3.6. Waves in the Isospeed Shallow 
Water Channel 


Basically, three methods were used to initially identify 
the arrivals on the digital oscilloscope. First, the 
arrival times predicted by the method of images were com- 
pared to the measured arrival times. The amplitude of 
the arrivals also helps in the identification process since 
the pressure amplitude is reduced by the reflection coeffi- 
cient at each of the interfaces of the layers. Finally, a 
path that reflects from the surface can be readily identi- 
fied by disturbing the water's surface with a ruler. Any 
surface-reflected arrivals change greatly on the oscillo- 
scope display and can easily be separated from those 


arrivals that only reflect from the bottom. 
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A typical display on the oscilloscope is shown in Fig. 
3.7. Here, the projector is inclined so that the beam has 
an incident angle of 55 dlagaes and the water depth is 35 cm. 
The arrivals shown are: 

1. direct wave, 

2. surface-reflected wave, 

3. bottom-reflected wave, 

4, wave reflected from the top of the plate, 


5. wave reflected from the bottom of the plate after 
passing through the plate, 


6. bottom-surface-reflected wave. 
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Figure 3.7. Different Kinds of Waves on Oscilloscope 
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Once the arrivals are identified for one source-to-receiver 
range, the signals can be followed by eye if the receiver 
Carriage is moved smoothly. 

Next, the range X is varied by moving the carriage and 
the arrival times T of the individual reflections are 
recorded. These measurements form the basis of the T-X and 
~?-x? methods. When these measurements are plotted as 
arrival time T as a function of range X, the direct arrivals 
appear as straight lines and the reflected arrivals (either 
from the top of the sand layer or from any of the layers 
below) appear as parabolas. (They are actually closer to 
hyperbolas but, over the limited range available in this 
experiment, parabolas were used as an approximation.) 

The T and X measurements were also plotted as Tt? as a 
hmoinection of xo On these coordinates, any of the lines are 
straight (or nearly so). From the slope of these lines, 
the sound speed either in the water or in the layers could 
be calculated. A linear regression method was used to 
compute the slope for the individual arrivals. 

Next, the ray parameter P was calculated. This value 
is fixed for a given measurement. The incident angle was 
determined from the source-receiver geometry and the T-X 
data were fit with parabolas. Then P could be calculated 
directly from the parabola equation: T = ax? + bX + Cc, 
aT/dX = P = 2aX + b. 

Once the parabolas have been determined, the layer depth 


can be found from the T-X measurements. For example, the 


zg 


arrival from the wave that passed through the upper sand 
layer would be described by a parabola of the form | 


Ty = a xt + bX, + cy) and the reflection from the top Of maa 
a xe + bee 


rae a, 2° 2: 


ray parameter P must be the same for each arrival; there- 


sand layer would be given by T The 


fore, for the wave passing through the sand layer, 


Pp = 2a,X) + by or Xy = (P-b,)/2a,, and, for the wave reflected 
from the top of the sand layer, P = 2a,X. + b, or X. = (P-b.)/2a.. 
2a 2 2 2 2 
Since AX — XK 7X5, | 
P=-b P=-b P=-b P-b : 
_ Li2 i - 22 B 
AT = [ia, EL AF by SE! AE cy} ta, ae rhe 


These two values (T and X) are necessary to calculate the 


depth and the sound speed (to compare with the t?-x? deter- 


mination of sound speed) of the layer. The values of layer 


depth Az and sound speed Cc, are given by, 


AEs 
ei yD. 
Az = 7 ¢,AT(L -P cy) 


The first value Of Cc.) Wseeconesene p?-x method and AT is 


1 
from the parabolic fit of the T-X data. The new value of 


c, can be compared to the old value as a check on the two 


al 


methods. 


C. SUMMARY OF EXPERIMENTAL PROCEDURE 
1. Bury the thin homogeneous plates at various depths; i 


2. Identify the arrivals on the oscilloscope; 


510 





Measure the horizontal distance X and the arrival 
time T of the various arrivals; 
Make the T-X and 7*-x? plots; 


Find the sound speed Cy in the layer by using the 
T“-x2 method; 


Determine the incident angle from the source- 
receiver geometry; 


Calculate the ray parameter P (fixed for each 
experiment) ; 


Fit the T-X data with parabolas; 
Calculate AX and AT by the ray parameter method; 
Calculate the layer depth Az; 


Calculate c, and compare to the value determined 
previously. 


sak 


IV. ANALYSIS AND RESULTS OF EXPERIMENTS 


A. MEASUREMENT AND ANALYSIS 


A measurement of the sound speed in the water was 


required for each of the reflection analyses so this was 


2 


measured using the T oe method. 


The slope of the straight 


line fit to the data points gave a sound speed of 1447 m/s 


and the correlation coefficient associated with the linear 


regression was almost 100%. 


For this measurement, the pro- 


jector beam was aimed horizontally and the hydrophone was 


positioned at the same depth as the projector. 


In each of the subsequent 
hydrophone depths were chosen 
and clarity of the reflection 
source Signal was four cycles 


1. Measurement 1 


experiments, the projector and 


to optimize the separation 


Signals. In every case, 


of a 150 kHz sine wave. 


the 


For this experiment an aluminum plate (thickness of 


1.7 cm) was buried 3.5 + 0.05 cm below the sand's surface 


and the projector depth was 8.5 cm while the hydrophone 


depth was 6.5 cm. 


TABLE I 


Speed of Sound (Measurement 1) 


WAVE SPEED OF SOUND (m/s) 
SAND 1661 
BOTTOM 1415 
Dike CE 1463 


32 


CORRELATION (3) 


ERROR (3%) 


Table I gives the speed of sound calculated from the various 
arrivals along with the correlation coefficients of the 
linear regressions. The speed of sound in the water should 
be compared to the value of 1447 m/s discussed above. The 
T-X and t-x? data are shown in Fig. 4.1. 

As discussed in the previous chapter, the T-X data 
were fit with a parabola for simplicity. The data for the 


wave that reflected from the top of the aluminum plate was 


fit with 


2 


Eoe= Om 2b xX — 2.89 X + ° 459 


and the data for the wave that reflected directly from the 


top of the sand was fit with 


T = 0.113 x* - 3.03 x + 410 
The geometry selected for the ray parameter determination 
gave an incident angle of 29 degrees which corresponds to a 


¢ s/m. This value of slope 


ray parameter P equal to 3.32 x10. 
locates a point on each of the two parabolas. MThe differ- 
ence in X values gives a AX of 5.47 cm. The parabolic form 
for the reflection from the top of the plate gives the time 
difference AT of 51.4 uS. 

Next, the layer thickness AZ was calculated from 


Equation (14). In this case, Az was determined to be 3.56 cm. 


This differs by only 0.3% from the measured value for the 
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overlying sand layer. (The burial depth of the plate gives 
the thickness of the top sand layer.) The value for the 
layer's sound speed was also calculated from the ray param- 
eter method but its value was not significantly different 
from that calculated by the Do=Ke method. 
2. Measurement 2 

For the second measurement, the hydrophone depth was 
6.4 cm, the projector depth was 8.5 cm and the aluminum plate 
(1.7 cm thick) was buried at a depth of 2.9 + 0.1 cm in the 
sand. The calculated sound speeds are given in Table II and 


the T-X and ee data are shown in Fig. 4.2. 


TABLE II 


Speed of Sound (Measurement 2) 


WAVE SPEED OF SOUND (m/s) CORRELATION (3%) ERROR (3%) 
SAND 1647 v4. 3 

BOTTOM 1419 Ne) sie io 
DELRECT 1437 Doi oal O69 


Figures 4.3 and 4.4 show representative displays from 
the oscilloscope at two of the ranges at which measurements 
were taken. 

The direct wave (A), the surface-reflected wave (B), 
the wave that reflects from the top of the sand (C), the 


wave that reflects from the top of the buried plate (D), and 
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Measurement 2 at X = 34.5 cm. 


Figure 4.4. Measurement 2 at X = 39.6 cm. 
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the wave that reflects from the top of the sand and aga 
from the water surface (E) are shown in these figures. The 
vertical line marks the beginning of the wave that is 
reflected from the top of the plate and the number gives the 
arrival time. 

In this case, the curve for the arrival reflected 
from the top of the plate is, 

T = 0.473 X* - 32.0 X + 1020 
and the curve for the arrival reflected from the top of the 


sand is, 


To =5s0 2136 x? =~ 5272 -X%)+ 47 


Here, the incident angle was 28.3 degrees so P was 3.30 x10 ° 


s/m. The ray parameter calculations gave 3.87 cm for AX and 
54 uS for dT so the layer depth Az was 3.73 cm. The actual 
sand layer thickness was 3.4 cm so the error was slightly 
less than 10%. The value for sound speed in the water by the 
ray parameter method was 1474 m/s which is about 103% high. 
3. Measurement 3 

In this measurement, an aluminum plate 0.1 cm thick 
was buried to a depth of 2.3 + 0.1 cm. The projectorsdeow 
was 8.5 cm and the hydrophone depth was 5.1 cm. The result- 

2 ee 


ing speed of sound measurements (determined by the T -X 


method) are summarized in Table III and the T-X data and the 
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T*-x? data are shown in Fig. 4.5. Figures 4.6 and 4.7 


show representative displays from the oscilloscope at two 


of the ranges at which measurements were taken. 


TABLE III 


Speed of Sound (Measurement 3) 


WAVE SPEED OF SOUND (m/s) CORRELATION (3) ERROR (3%) 
SAND 1649 95.2 

BOTTOM 1471 99.9 ie 7 
DIRECGI oO oi) ets 4.3 


The parabola for the arrivals reflecting from the 


top of the plate was determined to be 


T = 0.290 X* - 13.3 X + 580 
and the parabola for the arrivals directly reflected from 


the top of the sand is 


T = 0.0484 x* + 0.528 xX + 357 


For the optimum geoemtry, the incident angle was 27 degrees 
so that the ray parameter P was 3.00 xg 
The ray parameter analysis gave 2.53 cm for AX and 


33.38 us for AT; therefore, the sand layer depth was computed 


to be 2.42 cm: less than one percent from the actual value. 
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T-X and 7?.x? Plots (Measurement 3) 


Figure 4.5. 
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Figure 4.6. Measurement 3 at X = 21.7 cm. 


Figure 4.7. Measurement 3 at X = 29.3 cm. 


41 





The calculated value for the sound speed in the sand layer 
was 1578 m/s and this differs by only 3 percent from the 
value obtained by the T*-x? method. 
4. Measurement 4 

For this measurement, an aluminum plate 0.3 cm thick 
was buried to a depth of 1.8 + 0.1 cm. The projector depth 
was 8.5 cm and the hydrophone depth was 3.5 cm. The speed 
of sound measurements is summarized in Table IV and the T-X 
data and the ~*-x? data are shown in Fig. 4.8. Representative 


displays from the oscilloscope at two ranges are shown in 


Figures 4.9 and 4.10. 


TABLE LV 


Speed of Sound (Measurement 4) 


WAVE SPEED OF SOUND (m/s) CORRELATION (3%) ERROR (3%) 
SAND 1582 oo 

BOTTOM 1458 oo Op 
DIRECE 1461 99R ae 


The parabolic form for the reflections from yeheouee. 


of the pillatemis 


T = 0.0259 x* + 0.361 X + 429 


and the form for the reflections from the top of thevsancuare 
2 


T = 0.0328 X + J See 
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Figure 4.9. Measurement 4 at X = 35.1 cm. 


Figure 4.10. Measurement 4 at X = 47.7 cm. 
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For the optimum geometry, the incident angle was 50 degrees 
so that the ray parameter P was 5.3 ae. 

The ray parameter analysis gave 5.80 cm for AX and 
40.4 us for AT; therefore, the sand layer depth was computed 
to be 1.74 cm. The calculated value for the sound speed in 
the sand layer was 1647 m/s. In both cases, the results 
compare very well with the measured value for layer depth 
and the noe value for sound speed. 

5. Measurement 5 

For this measurement, the projector depth was 8.5 cm 
and the hydrophone depth was 3.5 cm. The aluminum plate 
(O.1 cm thick) was buried at a depth of 5.0 + 0.1cm. The 
speed of sound measurements are summarized in Table V and 
the T-X data and the p?_x? data are shown in Fig. 4.11. 


Samples of the received arrivals at two ranges are shown in 


Figures 4.12 and 4.13. 


TABLE V 


Speed of Sound (Measurement 5) 


WAVE SPEED OF SOUND (m/s) CORRELATION (3) ERROR (3%) 
SAND ey yy 9959 

BOTTOM 1468 100 1.4 
PERE CT 1488 100 2,8 


The parabolic form for the reflections from the top 


of the plate is 
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T = 0.0450 x* - 0.213 X + 407 
and the form for the reflections from the Gopmersenen-om 
is 


T = 0.0463 x* + 0.562 X + 339 


The incident angle for determining the ray parameter was 
33.4 degrees so that the ray parameter P was 3.70 rai O He 

The ray parameter analysis gave 9.59 emetore. 9a] 
71.5 us for AT; therefore, the sand layer depth iwaaweon= 
puted to be 4.58 cm--about six percent below the actual 
value. The calculated value for the sound speed in the sand 


layer was 1904 m/s which was significantly higher than any 


other of the determinations. 


B. OVERALL ORESULIS 

In all, seven experiments were performed but the data 
from two of the experiments was not good enough to analyZe. 
The other five experiments have been reported above and, 
except for number 5, the results were quite consistent with 
reasonable errors. 

Table VI gives a summary of the results of the five 
experiments that were analyzed and compares the results of 
the Toe method and the ray parameter method. The average 
error in determining the depth of the layer was 3.5 percent 


so this measurement can be done effectively; however, the 
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TABLE VI 


Summary of Brief Results 


MEASUREMENT 
CONTENT il 2 3 4 5 
4 ACTUAL Beor05)) Society 23 (e.1) 1,804.1) 5.0(+.1) 
4 EXPERIMENTAL 3.56 3.73 DEY, ee 4.58 
§  1°-x* METHOD 1662 1647 1649 1582 1577 
a RAY PARAMETER 1791 1474 1578 1647 1904 
ie DIFFERENCE 163 154 49.6 18.9 275.9 
© ERROR (8) 10.0 9.5 35.0 1.2 16.9 


error in the sound speed calculation (for the sand layer) 
was fairly large (about 8 percent). Consequently, the 
sound speed determination is not guite as reliable. 

In Table VI, the difference means the values between the 
average of the five p?-x? results and the ray parameter 


method's value. 
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V. DISCUSSIONS AND CONCLUSIONS 


The objective of this research was to localize buried 
objects (e.g., thin aluminum plates) in water-saturated sand. 
In order to accomplish this objective, we needed to predict 
the sound speed and the depth of the overlying layer. The 
results are summarized in Table VII. 

The average speed of sound by the Toe method is 1472 
m/sec for the direct wave, 1446 m/sec for the bottom reflected 
wave, 1623 m/sec for the wave through the sand. 

Comparing the above sound speeds with the value given by 
Reference 3 (1481 m/sec at 20 °C), we had 0.6% error in the 
direct wave speed and 2.3% error in the bottom reflected 
wave speed. 

Reference 3 gives the sound speeds 1540 m/sec for coarse 
silt and 1730 m/sec for quartz sand. The sand used in this 
experiment is #30 mesh sand (washed, kilndried and graded) 
and the measured sound speed is 1623 m/sec. 

We used linear regression for finding the speed of sound 


with the ~t_x2 method. Since most correlation values of the 


linear regression are over 95%, the Tepe method is reasonable 
for determining the speed of sound. 
By using the ray parameter method, we determined AX, AT 


and Az. The average error between the actual layer's depth 


and the computed depth is 3.5%. The error range is between 
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TABLE VII 


The Detailed Results of Total Experiment 


MEASUREMENT 1 2 S 4 > 


von) 3 OC eles (Sole (+.1)°95.0(+.1) 


98.3 94.3 95.2 99.7 99.9 
SPEED S 
‘i 1661 1647 1649 1582 1577 
re} iE 
O 
S SOUND 99.4 99.6 99.9 99.7 100 
= (m/sec) e 
~o 1416 1419 1471 1458 1468 
& 
j 
GN 
ae ; 99.8 99.9 99.8 99.9 100 
= CORRELA- 1463 1437 1510 1461 1488 
a GL Nt EE 
P (10° Sec /m) 3.32 3.69 3.00 5.30 3.70 
E AX (cm) 5.47 3.87 2.53 5.8 9.59 
ov 
‘ AT (usec) 51.4 54.0 33.9 40.4 71.5 
2 Az (cm) 3.56 3098 2.42 1.74 4.58 
8 ERROR (z)1 Cae 9.7 0.8 0.0 6.5 
g c, DIFFERENCE’ 163 154 50 19 276 
(m/sec) 
rb) 
A ERROR (3) HOC 9.5 3.0 2 16.9 


Average error Of measuring Az = 3.5% 


Average error of Cy difference = 8.1% 


a ae ep ee SS ee ee ee ee es es ee ee eS es ee es ee ee ee ee ee ee ee ee ec ee ee ES ee ee ee Se ee ee ee 


tihis error is the value when comparing the experimental Az with actual 
layer's depth. 


25 difference is the difference between Cy Oi T?_x2 and the value of 


the ray parameter method. 


3this error is the _value when comparing the Cy difference with the 
average value of T -X2, 
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0% and 9.7%. The ray parameter method is one of the most 
useful for finding an unknown layer's depth. But, the average 
Cnror OF Cy between the Doe Ke and the ray parameter methods 
is relatively high (8.1%), 

In every measurement the direct wave's speed of sound 
is always faster than that of the bottom reflected wave. 
These values should be the same because both are measures of 
the sound speed in water. This discrepancy is probably a 
result of the difficulty in determining the true horizontal 
distance between the acoustic center of the projector and 
the hydrophone for different angles of the projector. In 
some cases, the difference in C7) between the 7*-x? and the 
ray parameter methods is over 10%. This may result from the 
choice of curve-fit for the T-X data. We tried several ways 
to find the best curve fit, but the computed results were 
quite different from the expected ones except for those from 
parabolwe Tac: 


We had some limitations on the measurements because of 


the use of a small tank. In particular, there were inter- 
fering reflections. So, some of the signals arrive at the 
Same time. That made the analysis of the layer speed of 


sound and depth very difficult. Also, the pulse length was 
too long to be able to resolve the thin layer. Projeeeen 
ringing limits the pulse length so that it cannot be made 
shorter than a certain duration. Because of these problems, 
this experiment would be easier to do ina larger body of 


water. 


SY 


If we do this experiment in the ocean instead of small 
water tank, then identifying the wave signals is much 
easier because of less multipath interference. In the ocean, 
we can try ASPER (airgun sonobuoy precision echo recorder) 
developed by the Hawaii Institute of Geophysics as one of 


the recommended methods [Ref. 1]. 
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